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We present preliminary results on the natural line width, T(D* + ), 
in the transition D* (2010) + -> D°tt + . The D° is reconstructed in two 
decay channels: D° —> K~tt + and D° —> K~7r + 7r~7r + . We use Monte 
Carlo to simulate our reconstruction resolution in the difference of the 
reconstructed D*+ and reconstructed D° masses. We fit a relativistic 
P-wave Breit-Wigner form convolved with the measured resolution, to- 
gether with a background function, to the mass difference distribution. 
For the decay mode D° -> K~7r + we obtain T = 83.5 ± 1.7 ± 1.2 keV and 
Am = m D * + — m D o = 145425.5 ± 0.6 ± 2.6 keV, where the quoted errors 
are statistical and systematic, respectively. For the D° —> K~n Jr Ti~Ti Jr 
mode, we do not quote a width measurement at this time, but mea- 
sure Am = m^>*+ — m D o = 145426.5 ± 0.4 ± 2.6 keV. Additionally, we 
present preliminary results from the decay mode A c (2880) + —> A+7r + 7r" 
for the mass, width and quasi-two-body branching ratios of the A c (2880) + 
baryon. We measure m (A c (2880) + ) - m (A c (2286) + ) = 594.9 ± 0.4 ± 0.3 
MeV and V (A c (2880) + ) = 6.2 ± 1.4 ± 1.5 MeV, where the errors are 
statistical and systematic, respectively. 
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1 D* (2010) + line width 

The D* + line width, T (D* + ), provides an experimental check of spectroscopic models, 
and is also related to the strong coupling gjj* Dtt • In the heavy-quark limit this coupling 
can, in turn, be related to the universal coupling g pQ. There is no direct experimental 
window on the corresponding coupling in the B system, Qb^b-ki since there is no phase 
space for the decay £>* — >■ Bit. However, the D and B systems can be related through 
g, and this opens an avenue to the calculation of gB*B^ DO 

We study D*+ —> D°tt+ decays for two D° decay modes, and extract values of the 
width and Am, the mass difference between the reconstructed D* + and reconstructed 
D° [2]. The detector resolution for Am is approximately 300 keV, which is larger than 
the previously measured value of T (D* + ) which was only approximately 90 keV. 
Therefore, to extract T (D* + ) we fit the Am signal with a relativistic Breit-Wigner 
(RBW), which is described in the next section, convolved with a resolution function 
taken from a GEANT4 Monte Carlo simulation (MC) of the detector response. The key 
to this analysis is that, while the core resolution is larger than the extracted width, 
the RBW tails allow us to use information far from the central signal region to define 
the width. The RBW distribution, shown in Eq. ^ below, contains the natural 
decay width r . 

The best measurement of the width to date is T = 96 ± 4 ± 22 keV by the CLEO 
collaboration, where the errors are statistical and systematic respectively [3j. The 
largest source of systematic uncertainty for that analysis was the lack of statistical 
significance of the tests used to investigate systematic effects. In this analysis, we have 
much larger sample sizes for all modes investigated, which allows us to make tight 
selections and to investigate sources of systematic uncertainty with high precision. 

1.1 Data Selection and Fit Procedure 

The measurements are based on an integrated luminosity of approximately 480 fb _1 
recorded at or near the T (AS) resonance by the BABAR detector at the asymmetric 
energy PEP-II e + e~ collider. The BABAR detector is described in detail elsewhere [4]. 
We reconstruct the decay D* + — >> D°tt+ in two prompt charm channels for the D°. 
The pion from decay is described as the "slow pion" because of the limited phase 
space available for this decay. The selection criteria were chosen to give a very good 
signal to background ratio in order to increase the sensitivity to the long RBW tails; 
we have not optimized them for statistical significance. 

Our fitting procedure involves two steps. The first step is to use a triple Gaussian 
and a threshold function to model the resolution due to track reconstruction by 
fitting the Am distribution using truth-matched Monte Carlo (MC) events. This 
MC was generated with a width of 0.1 keV, so that all of the observed spread is due 
to reconstruction effects. The threshold function includes a non-Gaussian component 
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that describes events in which the tt s decayed to a /x in flight and we use hits from both 
the 7r and /i segments. Accounting for these non-Gaussian events greatly improves 
the overall fit at low and high Am values relative to the signal peak. 

We use a binned maximum likelihood method to fit the MC data with the function: 

f th T (Am;p, a) + f ± G (Am; fi u a ± ) + f 2 G (Am; /x 2 , a 2 ) + (1 - fi - / 2 ) G (Am; /x 3 , a 3 ) 

(1) 

where the G(x;/^,c^) are Gaussian functions and fth> fi, f2 are the fractions con- 
tributed by the threshold function, and the first and second Gaussian components, 
respectively. The threshold function is defined as 

T (Am; p, a) = Am u p e a ' u (2) 

where u = (Am/Am ) 2 — 1 and a is a slope parameter. 

Figure [I] shows the individual resolution function fits for the two D° decay chan- 
nels. Each plot shows the total resolution probability density function (PDF) in 
green, the individual Gaussian contributions in magenta, and the threshold function 
in red. The resolution functions should peak at the generated mass value; however, 
there is a slight bias toward high Am in the weighted mean of the three Gaussians 
(the threshold function accounts for less than 1% of the signal and is excluded from 
this calculation). We take this bias as an offset when extracting the position of the 
RBW pole from data. The offset is only a few keV in size: 4.3 keV and 2.8 for Kn 
and X37T, respectively. 




Figure 1: The fit to the MC resolution data sample for both D° decay modes. The 
interval size is 50 keV, and the high mass tails are dominated by low statistics. 
Relative fractions for each PDF component are shown in the legend (rounded to the 
nearest 0.01). 

The second step is to take these resolution function parameters, convolve the 
Gaussian component with a relativistic Breit-Wigner (RBW) of the form given in 
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Eq. ([3]), and fit the Am distribution in data. We fix the shape parameters and 
relative fractions from the resolution function in the fit to data. There is additional 
flexibility in the fit to data as we allow a common offset, 5, for the Gaussian mean 
values and a width scale factor, e, between MC and data. 

dT(m) _ (p\ 2i+l ( l + rV \ (m r ) 2 
dm \ Po ) V 1 + r V J ( m § - m 2 ) 2 + (m r Total (m)) 2 

The second factor is from the Dtt form factor, represented by a Blatt-Weisskopf 
form factor with radius parameter r j5j [6]. We use the value r = 1.6 GeV -1 from 
Ref. [7J. The specific value of the radius is unknown, but for the charm sector it 
is expected to be ~ 1 GeV -1 . We vary this value as part of our investigation of 
systematic uncertainties. Since the D* + is a vector particle we use t — 1. r Tota i(m) 
is defined by Eq. Q and, subsequently, by Eq. ([5]). 

r-ibtaiM ~ r D * D ^(m) (4) 



rw») = r (iLY' +1 (5) 

Here, as in most applications, we approximate Y^ ^\{m) ~ r^L^m), by ignoring the 
electromagnetic contribution from D* + —> D +r y. This approximation has a negligible 
effect on the extracted values. We normalize the resulting PDF such that, at the pole 
m = m , Eq. ^ is equal to 1.0. Overall scale factors (such as factors of 7r, etc.) are 
swept into the normalization. Figure [2] shows the fits to data for both channels. The 
normalized residuals show very good agreement between the data and our model. 



1.2 Systematic uncertainties 

The systematic errors are quantified independently for the D° —> Ktt and L* — » K3tt 
modes, and are summarized in Table [T] 

For our studies of V and the RBW pole as a function of pi a ^ azimuthal angle 0, 
and reconstructed D° mass, we split the data into disjoint subsets. These particular 
variables were chosen for careful study based in on dependency issues observed in 
previous BABAR analyses. We ascribe an error for the variation observed in each 
variable using a procedure similar to the scale factor method used by the PDG (see 
their discussion of unconstrained averaging [8j). The width value does not display 
an azimuthal dependence; however, the RBW pole does. We continue to investigate 
the origin of this dependence, and temporarily take half of the variation as the error. 
Looking in slices of the reconstructed D° mass and pi^ we estimate small systematic 
errors on Y and the RBW pole. 
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Figure 2: The fit to data for both D° decay modes (shown with logarithmic y axis 
scale). The total PDF is shown in green, the signal component in blue, and the 
background in red. The Am interval size is 50 keV. The normalized residuals, shown 
below each figure, are defined as (N observed - N predicted ) / \/N ohserved . 



Our choice of fit range in the Am distribution is arbitrary, and so we study the 
effect of systematically varying the endpoint. The choice of radius in the Blatt- 
Weisskopf form factor was tested by reducing the description of the D*+ to a point 
particle. We also estimate the systematic uncertainties which result from the mod- 
eling of the detector material and from the measurement of the magnetic field. The 
resolution function parameters extracted from the MC fit are also varied. Since we 
use a binned maximum likelihood fit, we also double the number of intervals used in 
our fit as a check of the sensitivity to the nominal interval size of 50 keV. 

For our validation fits, we generate signal MC samples with values for T and 
m (D* + ) — m (D°) close to the current world averages [8j. Each signal MC sample is 
paired with a background cocktail in the same ratio as in the corresponding fits to 
data. The validation data samples are approximately 5 times larger than the data. 
For each D° decay mode, the validation fits describe the data very well. In the fit 
to K3n, however, we do observe a bias in the extracted width. To be conservative 
we take the full difference between the fitted and generated value of the width, and 
assign a 5.7 keV error. At present, we are investigating the source of this bias, which 
we hope to minimize. 
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Source 


a sys (r) [keV] 


<j s?/s (RBWpole) [keV] 


Kw 


K3n 


Kw 


K3n 


Detector material and magnetic field 


0.3 


1.3 


0.7 


0.75 


Plat, dependence 


0.57 


0.19 


0.2 


0.2 


m (D® n ) dependence 


0.0 


1.12 


0.01 


0.0 


Azimuthal dependence 


0.45 


0.65 


2.5 


2.5 


Am fit range 


0.8 


0.3 


0.1 


0.0 


Resolution shape parameters 


0.41 


0.3 


0.17 


0.14 


Blatt-Weisskopf radius 


0.04 


0.03 


0.0 


0.0 


Am interval size 


0.1 


0.1 


0.1 


0.1 


MC validation 


0.0 


5.7 


0.0 


0.0 


Total 


1.2 


6.0 


2.6 


2.6 



Table 1: Summary of systematic errors. Total error is the sum in quadrature of each 
column. 

1.3 Results 

The preliminary results of the fits are 

T(Ktt) = 83.5 ± 1.7 ± 1.2 keV 
Am (Kir) = 145425.5 ± 0.6 ± 2.6 keV 
Am (K3n) = 145426.5 ± 0.4 ± 2.6 keV 

At this time, we do not report the width measurement for the K3tt mode pending 
investigation of the bias in the validation fit. The measurements of Am for the 
two modes are in excellent agreement, and have small systematic uncertainties. The 
uncertainties reported here are conservative overestimates and we anticipate that they 
may be reduced further. 

The D* + line width can be related to the effective coupling of a pseudoscalar 
to a quark. The total line width includes contributions from the D* + —> D°tt + , 
D* + Z) + 7T°, and D* + — > D + ^f decay channels. The strong couplings gD*D ir+ an d 
9d*d+it are related through isospin conservation. Using chiral perturbation theory 
and the heavy-quark limit on the c quark, this D* + strong coupling can be related to a 
universal coupling constant, g. Di Pierro and Eichten obtain predictions for the value 
of this universal coupling constant in a variety of decay modes in terms of a ratio, i?, 
which involves the width of the initial state [9j. The coupling constant should then 
take the same value for the three D decay channels listed in Table [2j which shows 
the values of the ratio R extracted from the model and the experimental values for 
r, as they were in 2001. The values of g obtained were in reasonable agreement, 
and indeed, at the time of publication, g was consistent among all of the modes in 
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state 


width (r) 


R = T/g 2 (model) 


g 


D* (2010) + 


96 ± 4 ± 22 keV 


143 keV 


0.82 ±0.09 


Di (2420) u 


18.91^ MeV 


16 MeV 


i.09i:i? 


D* (2460) u 


23 ± 5 MeV 


38 MeV 


0.77 ±0.08 



Table 2: Selected rows from Table 11 of Ref. [9j. State names correspond to the 
current PDG listings. The third column is the ratio, R ) extracted from the model in 
Ref. [9j. The values of g were obtained from the data available in 2001. 



state 


width (r) 


R = T/g 2 (model) 


g 


D* (2010) + 


83.3 ± 1.3 ± 1.4 keV 


143 keV 


0.76 ±0.01 


D l (2420) u 


31.4 ± 0.5 ± 1.3 MeV 


16 MeV 


1.40 ±0.03 


D* (2460) u 


50.5 ±0.6 ±0.7 MeV 


38 MeV 


1.15 ±0.01 



Table 3: Updated coupling constant calculations using the latest width measurements. 
The values of R are as in Table [2} There is significant disagreement among the new 
values. 



Ref. [9j. In 2010, BABAR published much more precise results on the D\ (2460) and 
D\ (2460)° [10]. Using those results, our preliminary T (D* + ) measurement, and the 
ratios from Table [2] we can calculate new values for the coupling constant. Table 
[3] shows the updated results. We estimate the error on the coupling constant value 
assuming a^/T <C 1. The updated widths result in large disagreements among the 
calculated values of the universal coupling constant. 

2 Measurements of A c (2880) + A+tt+tt- 

A large number of charm baryons have been observed, but certain parameter values 
and decay characteristics still are not well-known [BJ. In this analysis, we focus on the 
A c (2880) + baryon, measure its mass and width, and search for resonant structures 
in its decay to the A c (2286) + 7r + 7r~ final state. The search yields the first mea- 
surements of quasi- two-body branching ratios of the decays A c (2880) + — >> Sc 7r+ an d 
A c (2880) + S++7T- where we reconstruct S c (2455) A c (2286) tt and S c (2520) -> 
A c (2286) 7r. Isospin symmetry demands equal branching fractions into the neutral 
and doubly charged modes. 
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2.1 Data Selection and Fit Procedure 



The measurements are based on an integrated luminosity of approximately 316 fb _1 
recorded at or near the T (AS) resonance by the BABAR detector at the asymmetric- 
energy PEP-II e + e~ collider. The BABAR detector is described in detail elsewhere 
[4]. The selection criteria are optimized using the S/^S + B metric on a 21.8 fb _1 
sample of A+ —> pK~n+ decays. This sample is excluded from the final analysis. To 

reduce B backgrounds we choose x p > 0.4 where x p = p*/p^ nax = V* I \j{ s l^) m \ c - 
Here, is the momentum in the center-of-mass (CM) frame, and s is the CM energy 
squared. 




Figure 3: AM distribution for the region x p > 0.8. The fit result is shown by the 
blue curve. 

To extract the mass and width, we fit the AM (= m (A+7r + 7r _ ) — m(A+)) dis- 
tribution with a RBW convolved with a Gaussian resolution function and a poly- 
nomial background function. The Gaussian width is fixed to the detector resolu- 
tion determined from MC. Figure [3] shows the fit result for the region x p > 0.8, 
which has the highest signal-to-background ratio. The preliminary results are AM = 
594.9±0.4±0.3 MeV and T = 6.2±1.4±1.5 MeV where the errors are statistical and 
systematic, respectively. The systematic error is the sum of the uncertainties associ- 
ated with the solenoidal component of the magnetic field, the tracking material, the 
background model, the signal model, the AM fit range, the x p momentum selection, 
and the detector resolution. 

Figure [| shows a number of resonant substructures, which are described in the 
caption. There are no visible overlaps in the A c (2880) + signal region indicated by the 
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M (A+7T+7T-) - M (A+) + 2.286 (GeV) 

Figure 4: Mass spectrum of A+7r + 7r~ invariant mass. Vertical bands correspond to 
A c (2593) + and A c (2625) + are clearly seen in the data. There is an enhancement in 
the A c (2765) + region. The red vertical lines enclose the relevant A c (2880) + range. 
Horizontal bands correspond to S c (2455) ++ and S c (2520) ++ decays. The small Q 
value of the S c decays create a kinematic correlation of the A c and 7r, which creates 
the diagonal bands corresponding to S c (2455)° and S c (2520)° decays. 

red lines, and so we ignore interference effects when we calculate the related quasi- 
two-body branching ratios. We fit the A c n projections with a Gaussian for the S c 
signal and a polynomial background. Then we obtain the event yield by performing 
a background subtraction using the sidebands defined for each decay mode (10 — 20 
MeV above and below the signal region). The raw yield is adjusted by an efficiency 
factor determined from MC to produce the corrected yield. Branching ratios are 
calculated according to Eq. ([6]). The non-resonant (NR) yield is the sum of the 
resonant yields subtracted from the total yield. 

B (A c (2880) + -> S c tt) _ Y (A c (2880) + S c tt) e (A c (2880) + -> S+tt+tt-) 
B (A c (2880) + ^ A+7T+7T-) ~ Y ( A c (2880) + ^ E+tt+tt") X e (A c (2880) + ^ S c tt) 

(6) 

Assuming no interference effects, and that BF (E c —¥ A c 7r) = 100%, the prelimi- 
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nary branching ratio values obtained are: 

B (S c (2455)° 7T+) jB (all A+tt+tt- 
B (E c (2520)° 7T+) jB (all A+tt+tt- 
B (£ c (2455) ++ 7T") /£ (all A+tt+tt- 
B (S c (2520) ++ 7T") IB (all A+tt+tt- 
£ (NRA+7T+7T-) IB (all A+tt+tt- 



3 Summary 

Our preliminary results of the D* (2010) + width and RBW pole are 

T(Ktt) = 83.5 ± 1.7 ± 1.2 keV 
Am (AV) = 145425.5 ± 0.6 ± 2.6 keV 
Am (AT37r) = 145426.5 ± 0.4 ± 2.6 keV 

At this time, we do not report the width measurement for the AT37T mode pending 
investigation of the bias in the validation fit. The measurements of Am for the 
two modes are in excellent agreement, and have small systematic uncertainties. The 
uncertainties reported here are conservative overestimates and we anticipate they 
may be reduced further. Using our result and updated width measurements for a few 
other states, we find large disagreements among the calculated values of the universal 
coupling constant using the ratios provided by Ref. 

Using data on the decay A c (2880) + — ¥ A c 7r + 7r~, we find preliminary values for 
the mass difference AM = m (A c (2880) + ) - m (A c (2286) + ) = 594.9 ± 0.4 ± 0.3 
MeV and the width T = 6.2 ± 1.4 ± 1.5 MeV. These values are consistent with the 
previous BABAR results from an analysis of the D°p mode [11] . We also report the first 
measurements of the branching ratios to the quasi-two-body E c 7r and non-resonant 
A + 7r + 7r~ states. 
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